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FOREWORD 


A  computer  program  is  presented  which  is  capable  of  modelling  the  geometry 
and  performance  of  the  water  entry  cavity  produced  by  the  entry  of  axially 
symmetric  projectiles.  The  program  is  based  on  an  existing  theoretical  cavity 
flow  model  modified  to  account  for  atmospheric  surface  pressure.  The  computer 
model  is  shown  to  satisfactorily  predict  cavities  produced  by  vertical  and 
oblique  entries  (including  high-speed  entries)  from  normal  air  of  right  circular 
cylinders  and  truncated  ogives. 

The  work  reported  herein  was  supported  by  NAVSEA  Code  63R31.  The  author 
would  like  to  acknowledge  Dr.  Thomaa  Peirce  of  NAVSEA  for  his  continued  interest 
and  support  of  water-entry  research.  The  author  also  extends  thanks  to 
Dr.  John  Baldwin  and  Messrs.  Charles  W.  Smith,  Jacob  Berezow  and  Robert  Kavetsky 
for  their  guidance,  advice  and  support.  a 
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INTRODUCTION 


The  first  part  of  this  report  presents  the  final  version  ot  a  computer 
program  which  is  capable  of  modelling  the  behavior  of  the  water  entry  cavity 
produced  at  vertical  water  entry.  The  program  is  based  on  a  hydraulic  cavi  ty 
flow  model  proposed  some  years  ago  by  G.  Birkhoff  and  K.  Isaacs^  which  hes 
been  modified  to  Include  Che  effects  of  an  atmospheric  surface  pressure. 

The  second  part  of  this  report  presents  a  summary  > the  results  of  a 
theoretical  and  experimental  study  of  the  water  entry  cavity.  In  this  study, 
observed  cavities  are  compared  with  cavities  predicted  by  the  computer  program 
described  in  Part  1.  Various  water  entry  conditions,  covering  a  wide  range  of 
entry  velocities,  were  investigated. 


iBirkhoff,  G.  and  Isaacs,  K.,  "Transient  Cavities  in  Air-Water  Entry," 
NAVORD  1490  (1951). 
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PART  1  -  CAVITY  MODELLING  CODE 


BACKGROUND 


Birkhoff  and  Isaacs  developed  two  flow  models  -  one  which  applies  Co 
vertical  entry  and  the  other  to  oblique  entry.  The  discussion  here  will  be 
limited  to  the  vertical  entry  model  on  which  the  program  is  based.  The  original 
model  is  applicable  to  the  prediction  of  the  open  cavity  produced  by  the  water 
entry  of  axially  symmetric  missiles  from  either  a  vacuum  or  greatly  reduced 
pressure  atmosphere  (Abelson2  found  that  it  also  gives  good  shape  predictions 
of  cavities  produced  by  normal  air  entry  up  until  the  time  of  pullaway). 

Birkhoff's  vertical  entry  model  is  based  on  the  hypothesis  that  flow 
"streamlines  lie  on  concentric  spherical  surfaces  centered  at  the  point  of 
impact."  Accordingly,  the  liquid  medium  can  be  thought  of  as  being  subdivided 
into  infinitesimally  thin  concentric  spherical  shells  or  "laminas"  centered  at 
the  entry  point,  as  shown  in  Figure  1.  As  the  missile  passes  a  lamina  the  fluid 
within  is  displaced,  or  forced,  radially  outward  as  shown  in  Figure  2.  The 
inside  edges  of  the  open  laminas  together  form  the  cavity  wall  outline.  The 
shells  derive  their  initial  kinetic  energy  from  the  loss  of  kinetic  energy  of 
the  passing  missile,  resulting  from  the  cavity  drag.  The  cavity  drag 
coefficient  of  the  missile  is  assumed  to  remain  constant  resulting  in  an 
exponential  velocity  decay  with  depth 


v 

P 


r 

P 


(1) 


where  • 

VpQ  *  missile  entry  velocity 
rp  ■  missile  depth 
a  *  retardation  coefficient 


2  Abel  son,  11.  ,  "The  Behavior  of  the  Cavity  Formed  by  a  Projectile  Entering  the 
Water  Vertically,"  Ph.D.  Thesis,  Mechanical  Engineering,  Univ.  of  Md.,  1969. 
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Considering  a  typical  lamina  located  a  distance  r  lrom  the  entry  point , 
Figure  3,  only  the  location  of  the  inside  edge  of  the  displaced  shell,  bb 
defined  by  r  and  the  angle  6g,  is  required  in  order  to  determine  the  position 
of  the  cavity  boundary.  By  treating  the  shell  as  a  kind  of  "fluid  pendulum"  and 
utilizing  the  principle  of  conservation  of  mechanical  energy,  the  differential 
equation  relating  the  angle  0a  to  time  can  be  shown  to  be'->- 


+  J 

^  (co* 

CD 

I 

coa  0ja) 

1  \v 

r  J  sin^  6^  In  F 

(2  - 

coa  6a) 

(1  +  coa  0a) 

cos 

ea  u  - 

coa  9a) 

(2) 


where  the  expression  for  cos  ^  is* 


cos  9m 


1 

2 


■a  tap  Vp2  e-2ar 
Tt  r3  pw  g 


(3) 


The  angle  9,,,  is  the  maximum  value  of  the  displacement  8a  corresponding  to 
the  peak  of  the  lamina's  swing.  The  initial  condition  can  be  found  by 


R 

(o)  “  erctan 
flo  a  r 


(4) 


where  0a  ,  the  initial  shell  angle,  is  taken  by  Birkhoff3  to  be  the  angle 
subtended  at  the  point  of  entry  by  the  zone  of  separation  of  the  water  from  the 
missile  surface.  This  would  imply  that  Kp  should  be  taken  bb  the  radius  of 
the  missile  at  the  zone  of  separation.  For  certain  missile  configurations  it  is 
obvious  where  the  separation  occurs  and  thus  what  Kp  should  be  (for  example, 
for  a  right  circular  cylinder  only  the  nose  flat  is  wetted  hence  Rp  would  be 
taken  as  half  the  nose  flat  diameter).  For  certain  other  missile  configurations, 
on  the  other  hand,  the  location  of  separation  may  not  be  obvious  and  may  be 
difficult  to  determine.  In  such  cases  it  may  be  necessary  to  estimate  Rp  (for 
example  in  Abelson's  cavity  study3  which  involves  sphere  entries,  Kp  is 
taken  as  the  transverse  radius  of  the  missile  body  at  its  center  of  gravity, 
i.e.,  hall  the  diameter  for  a  sphere).* 


3See  footnote  1  on  page  7. 

3Abel  son,  H. ,  "Cavity  Shapes  at  Vertical  Water  Entry  -  A  Comparison  of 
Calculated  and  Observed  Shapes,"  NOLTR  67-31  (1967). 

*ln  reference  (3)  the  expression  for  cos  0m  is  in  error  in  that  the  first  term 
in  the  brackets  in  Equation  (3)  appears  as  a  "l"  instead  of  a  "2"  as  it  should 
be.  The  author,  apparently  through  an  oversight,  failed  to  account  for  the  "at 
rest"  potential  energy  of  the  shell  which,  because  of  the  reference  datum 
selected,  is  not  equal  to  zero.  If  this  at  rest  term  is  included  in  the 
reference  (3)  derivation  the  correct  expression  for  cos  0m,  with  the  "2"  in  the 
brackets  instead  of  the  "1",  is  obtained. 
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By  differentiating  (2)  with  respect  to  time  the  expression  for  the  angular 
acceleration  is  obtained: 


6  _  "is 

a  r 


1  +  (cos  8  -  cos  8  )  Z 

_ a _ m 

sin  0  In  00 

La  J 


where : 


2-4  cos  9. 


2  cos  6 


■  —  -a  — 

cos  0  (cosJ  0  -  2  c o s 0  -  cos  8  +  2 )  In  (xT  sin^  0 

a  a  a  a  a 

(2  -  cos  0.)  (1  +  cos  0  ) 

^  —  _ _  ,  &  a  _ 

cos  8a  (1  -  cos  8a) 

MODIFICATIONS  TO  ORIGINAL  FLOW  MODEL 

The  foregoing  equations,  which  describe  Birkhoff’s  original  model,  do  not 
account  for  pressure  either  in  the  cavity  or  at  the  water  surface  since  they 
were  assumed  to  be  zero.  However,  as  mentioned  earlier,  the  model  as  used  in 
the  program  has  been  altered  to  include  the  effects  of  surface  pressure  since 
this  modification  is  fairly  straightforward,  providing  the  surface  pressure  is 
regarded  as  constant.  Modifications,  on  the  other  hand,  to  include  the  effects 
of  pressure  in  the  cavity  on  cavity  behavior  are  not  a  simple  matter  at  all.^ 
Thus,  in  the  program  cavity  pressure  is  assumed  to  equal  zero. 


When  the  effect  of  surface  pressure,  P0,  is  included  it  is  shown  in 
Appendix  2  that  the  governing  differential  equation  for  the  displacement  angle 


8a  becomes; 


&a  “  + 


4 

rg  +  CPq/Pw) 

COB  0  -  cos  0 

a  ta 

2  .2  Q  ..  rvi 

where,  as  in  the  original  model,  the  initial  condition  is 


0,,  *  arc  tun 

u  o 


with  corresponding  relations  for  9a  and  cos  ^  as 

-2  [rg  +  (P  /p  )|  I-1  +  (cos  0  -  cos  9  )  (Z)' 

fl  -  - SL-JLl - * - - 


sin  9  In 
a 


^See  footnote  2  on  page  9. 


11 


NSWC  TR  81-59 


and 


cos  0„ 


am 

2  -  -  -g 

TT  r2  (r 


P  e 
o _ 


-2ar 


Pw  8  +  Po) 


(8) 


It  can  be  seen  that  by  setting  P0  equal  to  zero,  Equations  (6),  (7)  and 
(8)  reduce,  to  Equations  (2),  (3)  and  (5). 


Clearly  the  differential  equations  for  the  angle  0a 
and  the  modified  model  are  non-linear  and  must  be  solved 
using  numerical  integration  since  analytic  solutions  for 
yet  to  be  found.  The  program,  in  fact,  solves  Equations 
numerical  integration. 


for  both  the  original 
for  0a  versus  time 
these  equations  have 
(6) ,  (7)  and  (8)  using 


CAVITY  MODELLING 


It  should  be  remembered  that  these  lamina  equations  apply  to  a  single  lamina 
at  a  distance  r  from  the  entry  point.  Hence,  the  motion  of  a  lamina  at  any  depth 
can  be  found  independent  of  the  laminas  at  the  other  depths.  For  this  reason  the 
entire  cavity  may  be  modelled  satisfactorily  with  the  use  of  a  relatively  few 
number  of  laminas  located  at  selected  depths,  as  shown  in  Figure  4.  The  cavity 
outline  at  any  given  time  may  be  obtained  easily  by  plotting  the  positions  of 
the  edges  of  the  laminas,  then  connecting  the  appropriate  points  to  get  the 
cavity  outline. 


Another  less  accurate,  but  more  convenient  method  of  describing  cavity 
shape,  which  will  represent  the  cavity  for  all  times,  involves  plotting  the 
cavity  diameter  at  the  selected  depths  versus  time  after  entry.  The  result  is  a 
set  of  arc-shaped  curves  like  the  ones  shown  in  Figure  8.  Obviously,  from  the 
figure,  deeper  laminas  start  their  motion  at  later  times.  It  can  easily  be 
shown  that  when  cavity  drag  coefficient  is  constant  the  time  for  the  missile  to 
reach  a  lamina  at  depth  r  is: 


t  =  (ear  -  l)/otVp 
ro 

From  Figure  5,  the  diameter  of  the  circular  opening  in  the  lamina  at  a  given 
time  is  found  by  the  relation 

D  *  2r  ain  0a  (9) 

There  is  an  approximation  being  made  here  however.  Although  the  diametet 
computed  from  Equation  (9)  is  exact,  it  is  actually  the  diameter  of  the  cavity 
not  at  the  depth  r  but  at  the  more  shallow  depth  r(l  -  cos  0a)  since,  from 
Figure  5,  the  edge  of  the  lamina  rises  up  as  the  lamina  swings  out.  For  small 
enough  values  of  the  angle  0fl  this  should  present  no  problem.  However,  even 
for  fairly  large  values  of  0a,  the  approximation  should  still  be  reasonable 
since,  in  general,  the  cavity  wall  is  nearly  vertical  meaning  the  change  in 
diameter  with  depth,  at  a  given  time,  is  not  that  great.  A  final  consequence  of 
using  this  technique  is  that  at  a  given  depth  r,  the  cavity  diameter  will  never 
exceed  the  limiting  value  2  r. 
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The  program  is  designed  to  make  it  convenient  to  implement  either  of  the 
two  modelling  techniques  mentioned.  The.  output  data  are  in  the  form  of 
"Spherical  Lamina  Performance  Data"  tables  where  each  table  corresponds  to  a 
lamina  at  some  user  specified  depth  and  contains  information  concerning  depth  r, 
8a ,  mid  Cavity  Diameter  versus  absolute  time  after  entry  as  well  as  other 
auxilliary  information  such  as  current  model  depths  and  angular  velocities  and 
accelerations  associated  with  the  shell  displacement  angle  0a. 

MINIMUM  DEPTH  RESTRICTIONS.  The  selection  of  depths  at  which  to  model  the 
cavity  cannot  be  made  entirely  arbitrarily  since  the  flow  model  equations  break 
down  in  the  region  near  the  water  surface. 


For  any  given  entry  there  is  a  lamina,  located  at  a  critical  depth, 
designated  rer^,  which  acquires  just  enough  energy  from  the  passing  missile  to 
rise  up  to  where  9a  is  equal  to  90°  (i.e.,  entry  conditions  are  such  that 
for  this  particular  lamina  6^,  is  equal  to  90°).  Needless  to  say,  all 
laminas  at  depths  less  than  rcr.  acquire  sufficient  energy  to  swing  up  beyond 
the  point  0a  *  90°.  In  physical  terms  such  a  lamina  will  continue  to  swing 
upward  until  it  flows  into  itself  at  the  point  0  a  ■  90°  as  depicted  in 
Figure  6.  By  setting  cos  8  in  Equation  (6)  equal  to  zero  (corresponding  to 
8m  ■  90°)  and  replacing  r  with  rcr  ,  the  following  formula  is  obtained 
which  can  be  solved  by  tr ial-and-error  for  the  first  critical  depth  rcr  : 


0  = 


1 

2 


2 


a  m 


2  -2arc 


- £ 

tt  r* 


Crl  ^rcri 


1 

«  *  V 


Inspection  of  Equations  (6)  through  (8)  will  show  that  it  is  safe  to  model 
the  cavity  at  depths  greater  than  rcrp  which  assures  that  the  displacement 
angle  8a  will  not  exceed  90°.  Attempting  to  model  at  depths  equal  to  or  less 
than  rcri  will  result  in  a  singularity  in  Equations  (6)  through  (8)  when  0a 
reaches  90°,  although  modelling  up  until  this  condition  occurs  may  be  possible. 


A  second- critical  depth,  designated  rcv  ,  occurs  nearer  to  the  surface 
than  r cr ^ ,  and  is  the  depth  below  which  Equation  (8),  for  cos  0,^,  breaks 
down  completely.  Attempting  to  model  at  depths  less  than  rcr  will  result  in  a 
value  for  cos  0m  which  is  less  than  -1.  Even  so,  as  in  the  case  of  the  first 
critical  depth,  Equations  (6)  through  (8)  still  "work"  until  0a  reaches  90° 
(where  the  singularity  occurs)  although  it  is  not  known  at  this  time  if  the 
results  produced  will  be  meaningful. 


Figure  7  summarizes  the  minimum  depth  restrictions.  Since  these 
restrictions  apply  to  depths  in  a  Bmall  region  near  the  water  surface  they 
should  not,  in  general,  interfere  with  the  modelling  of  the  overall  cavity. 


PROGRAM  DESCRIPTION 


The  computer  program  is  written  in  standard  FORTRAN 
appears  in  Appendix  1.  The  first  part  of  the  listing  is 
manual."  To  demonstrate  the  use  of  the  program  a  sample 
carried  out  in  the  next  article.  A  brief  description  of 
here . 


and  a  listing  of  it 
essentially  a  "user's 
problem  has  been 
the  program  is  given 
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Basically  the  program  solves  the  lamina  motion  Equation  (6),  (7)  and  (8) 
for  0a  with  time  for  the  lamina  at  each  user  specified  depth  using  the  three 
term  taylor  series: 


0 


a( i  +  1) 


0“(i)  *  \i) 


*  2  8*(i,  <**>* 


Currently,  the  maximum  allowable  number  of  depths  which  may  be  specified  is 
10,  however  this  number  may  be  increased  simply  by  increasing  the  Fortran 
DIMENSION  parameters  for  the  variables  DEPTH,  DELTA  and  NSTORE  as  required.  The 
end  result  of  the  computations  for  a  lamina  at  a  particular  depth  is  the  single 
"Spherical  Lamina  Performance  Data"  table  described  in  the  previous  section. 
There  are  no  "Max  Time"  or  "Max  Iterations"  parameters  to  bother  with  since, 
when  the  motion  of  a  particular  lamina  is  completed,  the  sign  of  the  angle  0a 
goes  negative,  and  the  computations  for  that  lamina  are  terminated 
automatically.  A  sample  of  the  data  tables  appears  in  Appendix  1  following  the 
program  listing.  These  tables  are  the  results  for  the  sample  problem  described 
next . 


SAMI LE  PROBLEM 


The  program  was  used  to  model  the  cavity  produced  by  the  water  entry  of  a 
steel  right  circular  cylinder,  1.50"  diameter  by  2.15"  long,  from  normal  air,  at 
an  entry  speed  of  250  feet/second.  The  cavity  was  modelled  with  laminas  at 
depths  of  10,  20  and  30  calibers.  Several  test  runs  were  carried  out,  using 
progressively  smaller  values  for  the  time  step,  until  convergence  to  the 
solutions  for  all  three  laminas  was  obtained.  The  results  of  the  final  run  can 
be  found  in  Appendix  1  following  the  program  listing.  The  first  page  of  the 
output  is  a  Reference  Information  table  listing  data  describing  the  entry 
conditions.  Following  the  Reference  Information  are  the  three  "Spherical  Lamina 
Performance  Data"  tables  for  the  laminas  at  each  of  the  selected  depths.  The 
data  in  these  tables  were  used  to  generate  the  curves,  in  Figure  8,  of  cavity 
diameter  history  at  each  of  the  three  depths. 
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PART  2  -  CAVITY  STUDY 


BACKGROUND 

G.  birkhoff*  compared  cavities  predicted  by  hia  original  model  ior 
vertical  entry  with  experimental  data  concerning  a  one  inch  diameter  steel 
sphere  entering  water  vertically  from  a  reduced  pressure  atmosphere  at  a 
velocity  of  73  feet/second  with  a  surface  pressure  of  1/27  atmosphere. 

Abelson^  compared  the  cavities  predicted  by  birkhofi's  model  to  those 
generated  by  the  vertical  entry  of  spheres  into  water  from  reduced  pressure 
atmospheres.  The  spheres  had  diameters  of  either  0.5  or  1.0  inches  and  entry 
speeds  of  between  62  and  69  feet/sccond. 

Initially,  Birkhoff's  original  model  was  applied  only  to  open  cavities 
produced  by  entry  from  a  vacuum  or  reduced  pressure  atmosphere.  Then 
AbelBon^  compared  the  model  to  the  cavity  produced  by  a  normal  air  entry  and 
found  that  theoretical  and  observed  cavity  shapes  agreed  well  up  until  the  time 
of  pullaway.  Abelson  then  modified  the  model  to  include  the  effects  of  both 
constant  surface  pressure  and  a  measured  cavity  pressure.  He  then  compared  the 
cavities  predicted  by  the  altered  model  to  the  observed  normal  air  entry 
cavities  and  found  that  the  modified  model  predicted  cavity  shapes  which  were, 
in  fact,  considerably  smaller  in  size  than  both  the  observed  cavities  or  those 
predicted  by  the  original  model.  The  fact  that  Birkhoff's  original  model, 
unaltered  for  system  pressure,  satisfactorily  predicts  cavities  produced  by  boch 
normal  and  reduced  pressure  entries  suggested  to  Abelson  that  cavity  shape  might 
be  practically  independent  of  surface  or  cavity  pressure,  at  least  up  to  the 
time  of  pullaway.  In  an  attempt  to  confirm  this  hypothesis,  Abelson  compared 
the  cavities  produced  by  sphere  entries  from  both  normaL  and  reduced  pressure 
atmospheres  and  found  that  "the  cavities  were  identical  until  surface  closure 
and  then  began  to  deviate  increasingly  until  pullaway,  with  the  cavities  of  the 
reduced  pressure  entries  slightly  larger."  As  a  result  of  his  observations, 
Abelson  concluded  that  the  cavity  shape  is  in  fact  practically  independent  of 
surface  or  cavity  pressure  up  to  the  time  of  pullaway  and  hence,  that  Birkhoff's 
original  model  is  the  one  to  use  to  satisfactorily  predict  cavity  shape  and  that 
"any  modification  of  the  model  to  include  system  pressure  invalidates  it." 


^See  footnote  1  on  page  7. 

n 

See  footnote  2  on  page  9. 

2 

See  footnote  3  on  page  10. 
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EXPERIMENTAL/THEQRBTICAL  STUDY 

Comparison!  of  Birkhoff's  vertical  entry  model  to  experiment  are  limited 
mostly  to  shots  involving  spheres  entering  at  relatively  low  velocities  (below 
73  feet/second).  An  exception  is  Abelson's  comparison  of  the  model  to  the  case 
of  a  water  entry  involving  a  3  inch  diameter  cylindrical  missile  having  a  140 
degree  conical  nose,  entering  at  200  feet/aecond  from  normal  air. 2  Jn  view  of 
these  limitations  it  was  felt  that  it  would  be  useful  to  compare  both  the 
original  model  and  the  version  modified  for  surface  pressure  (in  spite  of 
Abelson's  conclusion  that  inclusion  of  system  pressure  invalidates  the  model)  to 
entries  involving  higher  velocities  and  model  configurations  other  than  spheres. 

Specifically,  theoretical  cavities  predicted  by  the  models  were  compared  to 
observed  cavities  from  a  number  of  shots  involving  normal  air  pressure  entries 
of  right  circular  cylinders  of  various  diameters  and  aspect  ratios  and  also  to 
shots  involving  cylindrical  missiles  of  various  sizes  having  truncated  ogive 
nose  configurations. 

One  source  of  experimental  data  was  photographic  records  from  vertical 
water  entry  teats  conducted  in  the  NSWC  Pilot  Tank  Facility  in  1967  under  the 
supervision  of  Dr.  A.  May.**  The  shots  chosen  for  this  Btudy  involved  right 
circular  cylinders  and  stepped  right  cylinders.*  Model  and  test  condition  data 
for  the  shots  chosen  from  May's  test  for  this  study  are  given  in  Table  1.  Nose 
flat  diameters  ranged  from  0.9  to  1.5  inches  and  entry  velocities  ranged  from 
175  to  415  feet/second. 

High-speed  motion  pictures  from  another  water  entry  test  program,  conducted 
in  the  NSWC  Undersea  Weapons  Tank  around  1959,  provided  a  second  source  of 
experimental  data.^  This  test  scries  involved  high  speed  oblique  entries  of 
truncated  ogives  and  right  circular  cylinders.  Although  these  shots  involved 
oblique  entries  it  was  felt  that,  for  the  depths  of  interest  in  this  study 
(below  50  calibers**),  the  entry  angles  of  45,  47  and  70  degrees  (from 
horizontal)  were  steep  enough  that  changes  in  the  cavity  geometry,  as  compared 
to  a  vertical  entry,  would  not  be  appreciable.  For  the  shots  chosen  from  the 
series  for  this  study  the  entry  model  nose  flat  diameters  ranged  from  0.707  to 
1.570  inches  and  entry  velocities  ranged  from  426  to  1175  feet/second.  Model 
and  test  condition  data  are  given  in  Table  2. 


^See  footnote  2  on  page  9. 

**May,  A.,  "The  Cavity  After  Vertical  Water  Entry,"  NOLTR  68-114  (1968). 

^May ,  A.  and  Hoover,  W.  R. ,  "A  Study, of  the  Water-Entry  Cavity," 

NOLTR  63-264  (1963). 

*For  the  two  shots  involving  the  stepped  cylinders,  the  stepped  shoulder 
remained  within  the  cavity  and  was  not  observed  to  interfare  in  any  way 
with  cavity  production. 

**Based  on  noae  flat  diameters. 
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TABLE  1 


MODEL  AND 

TEST  CONDITION  DATA: 

PILOT  TANK 

SERIES 

MODEL 

NOSE  FLAT 

MODEL 

MODEL 

SHOT  # 

CONFIG. 

DIA. 

LENGTH 

WT. 

MAT'L. 

(in)  ( in)  TlbeT 


1419 

RIGHT  CIRC.  CYL. 

1.000 

8.25 

1.834 

STEEL 

1711 

RIGHT  CIRC.  CYL. 

1.500 

2.00 

2.143 

UNKNOWN 

1438 

RIGHT  CIRC.  CYL. 

1.500 

4.39 

2.195 

STEEL 

1434 

RIGHT  CIRC.  CYL. 

1.500 

2.15 

1.075 

STEEL 

1474 

SINGLE  STEPPED  CYL.* 

0.900 

10.10** 

1.818 

STEEL 

1476 

SINGLE  STEPPED  CYL.* 

0.900 

10.10** 

1.818 

STEEL 

ENTRY 

DRAG 

CAVITY 

ENTRY 

SHOT  # 

VELOC. 

AREA 

DRAG  COEFF. 

ANGLE 

( ft/ sec)  (sq  in) 


1419 

175.0 

.785 

.807 

VERTICAL 

1711 

200.3 

1.767 

.807 

VERTICAL 

1438 

245.0 

1.767 

.807 

VERTICAL 

1434 

250.0 

1.767 

.807 

VERTICAL 

1474 

320.0 

0.636 

.807 

VERTICAL 

1476 

415.0 

0.636 

.807 

VERTICAL 

♦Stepped  cylinder  dimensions:  0.90"  dia  x  0.975"  long  nose  cylinder; 

1.50"  dia.  x  3.308"  long  body  cylinder. 

♦♦Length  for  an  equivalent  steel  unstepped  cylinder  (i.e.,  one  having  the 
same  weight  and  nose  flat  diameter  of  the  stepped  cylinder). 
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TABLE  2 


MODEL  AND  TEST 

CONDITION 

DATA:  UNDERSEA 

WEAPONS  TANK 

SERIES 

MODEL 

NOSE  FLAT 

MODEL 

MOD 

SHOT 

1 

CONFIG . 

DIA.* 

LENGTH 

WT. 

MAT 1 L . 

(in) 

lin) 

(lbs) 

184 

TRUNCATED  OGIVE 

0.707" 

NOT  AVAILABLE 

1.599 

NOT  AVAILABLE 

192 

TRUNCATED  OGIVE 

0.864 

NOT  AVAILABLE 

1.600 

NOT  AVAILABLE 

205 

TRUNCATED  OGIVE 

.774 

NOT  AVAILABLE 

1.669 

NOT  AVAILABLE 

193 

TRUNCATED  OGIVE 

.864 

NOT  AVAILABLE 

1.600 

NOT  AVAILABLE 

63 

TRUNCATED  OGIVE 

.785 

NOT  AVAILABLE 

1.652 

NOT  AVAILABLE 

172 

RIGHT  CIRC.  CYL. 

1.570 

NOT  AVAILABLE 

1.537 

NOT  AVAILABLE 

ENTRY 

DRAG 

CAVITY 

ENTRY 

SHOT 

1 

VELOC. 

AREA 

DRAG  COEFF. 

ANGLE** 

, 

(f t/aec) 

(sq  in) 

(deg) 

184 

426.0 

.393 

.800 

47.0 

192 

548.0 

.586 

.800 

47.0 

205 

599.0 

.471 

.800 

47.0 

193 

832.0 

.586 

.800 

47.0 

63 

946.0 

.484 

.800 

70.0 

172 

1175.0 

1.936 

.800 

45.0 

*For  these  models  full  body  diameters  were  equal  to  1.57  inches. 

**Angle  between  theoretical  gun  line  and  water  surface. 
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SUMMARY  OF  RESULTS 

All  together,  tweLve  water  entries  were  modelled  with  results  presented  in 
the  form  of  cavity  diameter  histories  shown  in  Figures  9  through  23.  The  curves 
in  Figures  9  through  20  compare  observed  cavities  with  cavities  predicted  by  the 
altered  theoretical  model  which  includes  effects  of  surface  pressure.  The 
remaining  curves  in  Figures  21  through  23  compare  observed  cavities  with 
cavities  predicted  by  Birkhoff's  original  model.  Each  graph  contains  cavity 
diameter  time  histories  for  both  predicted  and  observed  cavities  at  selected 
depths.  The  data  in  nearly  all  the  plots  extend  well  beyond  the  time  of 
pullaway  (the  only  exception  is  shot  #172  where  a  deep  closure  occurred  and  no 
pullaway  was  observed).  Time  data  were  obtained  from  timing  marks  present  on 
th  e  film. 

For  the  curves  in  Figures  9  through  23,  corresponding  to  the  study  with 
surface  pressure,  the  following  observations  were  made: 

1.  In  general  the  model  altered  for  surface  pressure  appears  to 
underestimate  the  cavity  diameter,  hence  cavity  size.  This  observation  should 
be  expected  to  some  degree  since  cavity  presure  is  assumed  to  be  zero.  The 
exceptions  are  shots  #63  and  #1476.  In  both  these  cases  eurveB  for  observed  and 
predicted  cavities  agree  closely  for  all  depths  for  roughly  the  first  half  of 
the  curves,  corresponding  to  cavity  growth,  then  for  the  laBt  half, 
corresponding  to  cavity  collapse,  the  model  overestimates  the  diameters. 

2.  Model  correctly  predicts  a  deep  closure  type  of  collapse  (collapse 
from  the  bottom  up)  for  shot  #172  and  base  closure  type  collapse  (collapse  from 
top  of  the  cavity  downward)  for  the  remaining  shots. 

3.  Theoretical  "lag  times"  (times  needed  for  missile  to  reach  selected 
depths)  which  have  nothing  to  do  with  the  model  but  result  from  the  assumption 
of  constant  cavity  drag  coefficient,  agree  well  with  observed  lag  times. 

Following  the  cavity  diameter  histories  is  a  set  of  curves  in  Figure  24 
showing  the  results  of  a  parametric  study  of  the  surface  pressure  P0  to 
determine  the  sensitivity  of  cavity  size,  at  a  representative  depth  of  20 
calibers,  to  variations  in  surface  presoure. 

The  cavity  diameter  histories  in  Figure  25  are  for  a  theoretical  cavity 
predicted  by  Birkhoff's  original  model.  It  is  included  to  show  how  the  typical 
theoretical  cavity  behaves  since,  for  the  study  involving  the  original  model  in 
Figures  21  through  23  only  a  small  portion  of  the  theoretical  curves  could  be 
shown  on  the  graphs  because  of  their  extremely  long  periods,  as  compared  to  the 
experimental  curves.  The  diameter  histories  in  Figure  25  (essentially  an 
extension  of  the  theoretical  curves  in  Figure  15  for  shot  #1474)  indicate  that  a 
deep  closure  occurs  at  a  depth  of  140  calibers. 

RECOMMENDATIONS 

The  following  recommendations  suggest  areas  to  consider  for  future  work  and 
items  which  may  be  of  concern  in  a  more  involved  study: 
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1.  Some  years  ago,  according  Co  H.  K.  Steves,  there  were  ahota  done  in 
the  NSWC  Hydro-Balliatica  Tank  Facility  involving  high  speed  vertical  entries 
which  are  recorded  on  film.  It  may  be  useful  to  see  how  well  cavities  produced 
by  the  model  compare  with  these  cavities  since  the  high  speed  shots  considered 
in  this  study  involved  oblique  entries  whereas  the  flow  model  is  intended  to  be 
applied  to  vertical  entries. 

2.  Reread  the  photographic  records  for  the  vertical  entries  in  Figures 
21  through  23  to  determine  approximately  when  the  time  of  pullaway  occurs  and 
indicate  these  times  on  the  graphs.  Then  compare  the  results  in  these  curves 
with  Abelson's  suggestion  that  cavity  shape  can  be  satisfactorily  predicted  up 
until  pullaway  for  normal  air  entries  using  Birkhoff's  original  model. 

3.  Modif icati a  u  nould  ba  made  to  the  current  version  of  the  cavity 
modelling  code  to  enable  it  to  handle  the  problems  associated  with  modelling  the 
cavity  in  the  critical  minimum  depth  region.  Currently  no  such  capability 
exists  in  the  code,  consequently  programming  difficulties  related  to  the 
numerical  complications  described  earlier  will  arise  if  an  attempt  is  made  to 
model  in  this  region.  At  the  very  least,  changes  could  be  made  to  the  code  so 
that  it  will  simply  ignore  specified  depths  which  fall  within  the  critical 
region  and,  in  addition,  issue  a  message  warning  the  user  of  the  condition. 

More  sophisticated  modifications  could  allow  some  limited  modelling  within  this 
critical  region. 
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LINE 


FIGURE  1  LIQUID  MEDIUM  SHOWN  SUB-DIVIDED  INTO 
CONCENTRIC  SPHERICAL  SHELL  OR 
"LAMINA"  REGIONS 


L.O.F. 


FIGURE  2  DISPLACEMENT  OF  FLUID  SHELLS  TO  FORM  CAVITY 
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LAMINA  I 


FIGURE  4  MODELLING  THE  CAVITY  WITH  THE  USE  OF  LAMINAS  AT 
SELECTED  DEPTHS 


FIGURE  5  CAVITY  DIAMETER  APPROXIMATION.  DIAMETER  AT  DEPTH  r 
APPROXIMATED  BY  PROJECTING  D  DOWN  TO  DEPTH  r 


FIGURE  8  FLUID  SHELL  IMPACTING  ITSELF 


N8WC  TR  81-B9 


L.O.F. 


FIGURE  7  SUMMARY  OF  CRITICAL  MINIMUM  DEPTH  RESTRICTIONS 
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SAMPLE  RUN 

6NTVEL-  260  fpi 
GUN  ANG  -  VERTICAL 

MODEL  GEOM  -  1.6"  OIA  x  2.16"  LG  RT  CYL  (STEEL) 
MODEL  WT- 0.904  LBS 
1  CALIBER  -  1  NOSE  FLAT  DIA 


FIGURE  B  CAVITY  DIAMETER  AT  VARIOUS  DEPTHS  VS  TIME 
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FIGURE  9  CAVITY  DIAMETER  AT  VARIOUS  DEPTHS  VS  TIME 


NSWC  TR  SI— 89 


NSWC  TR  81-50 


SHOT  #1711 

MODEL  OEOM  -  RIGHT  CIRC  CVL 
1  CALIBER  -  1  NOSE  FLT  DIA 
P0  -  1  ATM 

- — -  -  EXPERIMENTAL 

- —  -  THEORETICAL 


FIGURE  16  CAVITY  DIAMETER  AT  VARIOUS  DEPTHS  VS  TIME 
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SHOT  #1438 

MODEL  QEOM  -  RIGHT  CIRC  CYL 
1  CALIBER  -  1  NOSE  FLAT  DIA 
P„  -  1  ATM 

- -  EXPERIMENTAL 


FIGURE  17  CAVITY  DIAMETER  AT  VARIOUS  DEPTHS  VS  TIME 
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SHOT  #1434 

MODEL  QEOM  -  RIGHT  CIRC  CYL 
1  CALIBER  -  1  NOSE  FLAT  DIA 
P0  -  1  ATM 

- -EXPERIMENTAL 

- -  THEORETICAL 


0.010  0.020  0.030  0.040 


TIME  (SEC) 


FIGURE  18  CAVITY  DIAMETER  AT  VARIOUS  DEPTHS  VS  TIME 


SHOT  #1474 

MODEL GEOM  =  STEPPED  CYL 
1  CALIBER  =  1  NOSE  FLAT  DIA 


FIGURE  19  CAVITY  DIAMETER  AT  VARIOUS  DEPTHS  VS  TIME 
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(suaanvo)  ua.iawvia  axiavd 


FIGURE  20  CAVITY  DIAMETER  AT  VARIOUS  DEPTHS  VS  TIME 


#1419 


FIGURE  22  CAVITY  DIAMETER  AT  VARIOUS  DEPTHS  VS  TIME 


TIME  (SEC)  r 

DIAMETER  AT  VARIOUS  DEPTHS  VS  TIME 


=  0  lbs/>q  ft 


FIGURE  24  CAVITY  DIAMETER  AT  20  CAL  DEPTH  VS  TIME.  FOR  VARIOUS  SURFACE  PRESSURES 
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TERMS 


Missile  cross-sectional  drag  area 
Lamina  rim  surface  area 

Cavity  drag  coefficient  of  projectile  based  on  Ac 

Diameter  of  lamina  opening;  diameter  of  cavity  at  a  given  depth; 
missile  nose-flat  diameter 

Gravitational  acceleration 

Projectile  mass 

Atmospheric  surface  pressure 

Radius  (and  "depth")  of  spherical  lamina 

First  critical  minimum  depth 

Second  critical  minimum  depth 

Projectile  depth 

Diameter  of  axis-symmetric  projectile  body  at  zone  of  separation 

Elapsed  time  for  miBsile  to  reach  depth  r,  measured  from  entry 

Kinetic  energy  of  fluid  in  shell 

Kinetic  energy  of  projectile 

Projectile  velocity 

Projectile  velocity  at  water  impact 

Total  potential  energy  of  lamina/surface-pressure  system 
Potential  energy  of  lamina  due  to  gravity 


NSWC  TR  81-59 

TERMS  (Cont'd) 


V £n  Initial  potential  energy  of  lamina/surf ace-pressure  ayatem  in  the 

at-reat  (i.e.,  8a  ■  0)  poaition 

^peak  Potential  energy  of  lamina/aurface-presaure  ayatem  when  lamina  ie  at 

peak  of  it»  awing  (i.e.,  0a  -  9m) 

VPo  Potential  reaulting  from  atmoapheric  surface  pressure 

Wp  Work  done  by  pressure  P0  acting  on  lamina 

o 

X,  Z  Space  saving  substitution  parameters 

01  Retardation  coefficient,  a  “  C[>AcPw/2mp 

®a»  Fluid  shell  displacement  angles 

Initial  shell  angle 
o 

0ra  Maximum  possible  value  of  the  shell  angle  9a 

Pw  Density  of  water 
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73/T4  OPT* l 


FTN  4.6*452 


1 


5 


10 


15 


20 


25 


30 


35 


40 


45 


50 


55 


C 

C 

C 

C 

C 

c: 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

0 

c 

c 

c 

c 

c 

0 

c 

c 

c 

c 

c 

c 

c 

c 

c 

C 

L 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


-PRQGkAm  description 


THE  PROGRAM  LISTED  BELOW  IS  DESIGNED  TO  CARRY  OUT  A 
NUMERICAL  SOLUTION  Of  THE  THEORETICAL  EQUATIONS  OF  MOTION 
WHICH  APPROXIMATELY  PREDICT  BEHAVIOR  OF  WATER  ENTRY  CAV¬ 
ITIES  FOR  VERTICAL  WATER  ENTRY. 

THE  EQUATIONS  ARE  TAKEN  FROM  A  HYDRAULIC  CAVITY 
FLOW  MOOEL.PPOPOStU  BY  G.  B IHKHOFF  AND  R.  I SAACS ( "TRAN¬ 
SIENT  CAVITIES  IN  AIk  WATER  ENTRY,"  NAVOHO  1490  ,  1951), 
WHICH  IS  BASED  ON  THE  APPROXIMATE  HYPOTHESIS  THAT  STREAM¬ 
LINES  LIE  ON  CONCENTRIC  SPHERICAL  SURFACES  (OR  "LAM  I  NAS" ) 
CENTERED  AT  THE  POINT  OF  IMPACT. ANY  NOSE  CONFIGURATION  FOR 
WHICH  THE  DRAG  COtFF.  IS  KNOWN  MAY  Bfc.  MOOFLLED.THF  EQUA¬ 
TIONS  AS  USED  IN  THE  PROGRAM  HAVE  BEEN  MODIFIED  TO  InCLUOE 
AN  ATMOSPHERIC  PRESSURE  TERM. 

SPECIFICALLY  THE  PROGRAM  GENERATES  "l AMINA  PERFOR¬ 
MANCE  DATA  TABLES"  FOR  THE  LAMINAS  AT  USER  SPECIFIED  DEPTHS 
(MEASUHED  ALONG  ThE  L INE-OF-F I  RE ) , BY  MtANS  OF  A  STEP-BY- 
STEP  INTEGRATION  OF  THE  LAMINA  MOTION  EONS . COMBINING  THE 
UAT A  IN  THE  TABLES  FOR  LAMINAS  AT  SEVERAL  DEPTHS  MAKES  IT 
POSSIBLE.  TU  OBTAIN  THE  PERFORMANCE  HISTORY  OF  THE  OVERALL 
CAVITY. THE  OUTPUT  UATA  ARE  DISPLAYED  IN  TERMS  OF  POLAR 
COORDINATE  HOWEVER  FOR  CONVENIENCE  AN  APPROXIMATE  CAVITY 
UIAMfcTER  is  COMPUlEU  ANO  DISPLAYED  BUT  IS  ACCURATE  ONLY 
FOR  SUFFICIENTLY  SMALL  VALUES  OF  THE  VARIABLE  "ANGLE". 

DIRECTIONS  TU  SET  UP  THE  INPUT  DECK  APPEAR  BELOW 
ALONG  WITH  A  DESCRIPTION  OF  IMPORTANT  V AR I ABLES . THE  PRO¬ 
GRAM  FEATURES  A  PRINT  FREQUENCY  PARAMETER . "NSTORE" . WH I CH 
AIDS  IN  RtDUCING  I  ME  QUANTITY  OF  OUTPUT  PRODUCED . ALSU 
THE  USER  NEED  NOT  BOTHER  WITH  "MAX  ITERATION"  OR  "MAX  TIME" 

parameters  sjnce. when  the  motion  of  a  particular  lamina 

IS  COMPLETED. COMPUTATIONS  FOR  HAT  LAMINA  ARE  TERMINATED 
AUTOMATICALLY. 

PROGRAM  AUTHUR-  M.  A.  MET/GER  ,  SEPT.  1  9Hn 


*  INPUT  DATA  4 


FIRST  CARD  - 

VARIABLES  *  REFOIA.RL NOTH, WEIGHT , RMORHO .DGAREA.UGCOEF 
DESCRIPTION  *  REFERENCE  D l A M . ( I N , ) , MODEL  LE NOTH ( IN , ) , 

MODEL  WEIGHMlBS.)  ,MOD.  DENS.  (LB/CU.IN.  )  » 
DRAG  AREAISU.IN.) .DRAG  COEFF..RESP. 

FORMAT  ■  6F  1 0  •  5 


SECOND  CARD 


VARIABLES 

DESCRIPTION 

FORMAT 


■  tNT VEL » NOEP  TH 

»  LNTRY  VELOCITY (FPS) ,  NUMHER  OF  SPECI¬ 
FIED  DEPTHS  PARAMETER,  RESPECTIVELY. 

■  F 10.5*15 


A-2 
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FIN  4.6*452 


60 


65 


7  0 


76 


80 


86 


90 


C 

C 

L 

L 

L 

C 

c 

L 

C 

L 

L 

L 

c 

L 

L 

c 

L 

c 

c 

c 

l 

c 

L 

C 

L 

c 

c 

c 

l, 

L 

c 

L 

C 

c 

L 

c. 


TH 1 KO  St  I  - 

VARIABLES  «  UEPTH.DELTAiNSTORE 
DESCRIPTION  «  EACH  CARD  IV  THIS  SET  CONTAINS  A  USER 
SPECIFIED  DEPTH  ALONG  WITH  VALUES  FOR 
THE  DESIRED  TIME  INCREMENT  AND  PRINT 
t'KEUUENCY  FOR  THAT  DEPTH. DEPTHS  SHOULD 
bE  IN  UNITS  OF  CALIBERS  RASED  ON  REF 01  A. 
SUBMIT  ONE  OF  THESE  CARDS  PER  SPECIFIED 
DEPTH. 

FORMAT  “  2F  1 0  •  5  *  1 5 


*  VARIABLES  • 


VARIABLES  *  ANGLE ♦ANGVEL.ANOACC 
DESCRIPTION  ■  ANGLE  REPRESENTS  THE  ANGLE  SUHTENOtD  AT 
THE  POINT  OF  ENTRY  BY  THF  L INK-OF-F IRE 
AND  THE  CAVITY  WALL  AT  DEPTH,  ANGVEL 
AND  ANG ACC  ARE  ANGULAR  VELOCITY  ANU 

acceleration  associatted  with  ANGLE. 
TYPES  ■  REAL 


VARIABLES  *  WATRHO.PRESSH 

DESCRIPTION  »  WATER  DENSITY  AND  ATMOSPHERIC  PRf.SSUHE 
TERM  RESP. 

IYPES  =  HEAL 


95  L 

C 


•I*6*:-1 
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PROGRAM  CAVITY  73/74  OPT-1 


FTN  4,6+452 


PROGRAM  CAVITY  (INPUT, OUTPUT , TAPE5-I NPUT » TAPtb-OUTPUT ) 

C 

C 

C  *##•,«#**♦«»«##***#*♦#,#«,*»§***•»#**»#»*•#****»«#* 

c  preliminary  operations 

c  . . . . 

c 

REAL  MOMASS  .  MDLUfcP 

DIMENSION  DEPTH <10 ) , DEL T A  < 1 0 ) ,NSTORE(lO) 

GRAVTY-32.174 
WATRHO-1. 93945 
PRESSR-21 16«H 

HEAD (5* 1000)  HEFD1 A*HLN3TH*WE lUHT , RMDBHO » DO ARE A , DGCOEF 
READ (5,1010)  ENTVEL,NDEPTH 
DO  1  I*1,N0EPTH 

READ <5, 1020)  DEPTH < I ) ,DELT A ( I >  » NSTORE < I > 

1  CONTINUE 

C 

C 

C  *  PRINT  reference  info  taule  » 

C 

WRITE (6,2000)  REFDIA,RLNGTH»OGCOEF,DGAKEA.WEIGHT, 

.  RMOHHO,WATRHO»PRESSR«8RAVTY  ,E'NTVFL*NDEPTH 

REFDl A-REFOI A/12,0 
OGAHEA-OGAHEA/144,0 
MDMA9S-WEIDHT/32.174 

RETARD- <DGLOEF*DGARLA*WATRH0) / < 2 . 0*MOMASS ) 

WRIT*  (6*2005)  RETARD 

WRITE (6,20 10)  (DEPTH ( I ) , I »1 *NUEPTH) 

DCONV- 180,0/3, 141592 

C 

C 

C  »  PRINT  TABLt  HEAOIN3  * 

C 

WRITf (6,2020) 

L 

C 

C  **####*##4  »***#**)(,*##*##**####»,*  »###*«###*,«##**#* 

L  GENERATION  OF  LAMINA 

C  PERFORMANCE  DATA  TABLES 

C  *********** **##*«*# *#♦***«#**#** **##*#* 

L 

DO  50  I«1,NDFPTH 
C 
C 

C  *  INI  I'lALIZE  PARAMETERS  * 

C 

A TEMP-RE  T  ARD*MUMASS*ENTVEL*ENT VEL*EXP  ( -  ?•  0*RFTARD*DEE,TH  ( I )  * 
.  REFDIA) 

BTEMP»2.0*3,141S92*(DEPTH(I) ##2)#(REFD1 A»*2) 

CTEMP-DEPTH ( I ) »KEFDI A*WATRHO*GRAVT Y+PRFSSR 
COSMAX-l . 0~A  TEMP/ <BTEMP*CTEM°) 

DEGANG-4C0S (COSMAX ) *180,0/3, 141592 

AHTIME- (ExP(RETAHO*I)EPTH( 1 ) *REFDI A  > -1 »  0 ) / <  RFT  ARD*ENT VEL  > 
WRITE  (6,2(130  I  DEPTH  ( I )  ,  DELTA  <  I  >  ,  NS  TORE  (  T  )  ,  DEGANG 
L 
C 


ao/i' 


A4 
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PHOoR AM  CAVITY  73/74  OPT*l 


KTN  4,6*452 


C 

C 


c 

c 

c 

c 


c 

c 

c 

c 


c 


c 

c 

c 

c 


c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 


t: 

e 

c 

c 


c 

c 

c 

c 


H  u 


*  COMMUTE  INITIAL  AN3LE  * 
anole«atam(i.o/(2.o*oe°th( n » > 

C0SAN6*C0S(ANGLli) 

SI NANG*5 IN ( ANGLfc ) 


“  COMMUTE  INllIAL  AN3.  V£L.  * 

Call  vh.oc  (pressr*watrho,gra\/ty,co$max.  anglf . 

OEP1HII)  , HE  ED  I  A  *  ANGVEL I 


*  COMMUTE  INllIAL  AN3.  ACC.  * 

Call  ACCEL <PRE55H,WATRNO*GHAVTY*COSMAX.ANGLF  * 

OEPTH(I) .REEDIA.AMGACC) 

MDLUEP® «LOG  <  HE  I  ARD*ENTVEL*ABT IMF* I  •  0 ) / (BET  AHD*REFOI A ) 
C  iVUIA*t>.n*DEPTH(l)  *S  l  N  (ANGLE  ) 


4  RADIAN  TO  DEGREE  CONVERSION  * 

DK 0 ANG= ANGLE *DCUNV 
UFGVEL  ANr,VEL«UCUNV 
DEGACC«ANGACC*iJCUNV 


*  PRINT  NEXT  LINF  OF  TABLE.  * 

WRITE.  (h,2040)  AbTlMEiCAVOlA, DEGANG, 

UtGVEL»l)EGACC,  MULUEP 

STEP  BY  slEP  INTEGRATION  OF 
MOTION  EONS  FOR  lamina  AT 
CURRENT  UEPTH(I) 

«»•*•**#**»•**•**#»*„*»»*«**•«*«»«•<»#«*«*«•»<>•«»«« 

OLoANGsaNGLE 
I  STORES') 

I S  TORE* I  STOWE* 1 


«  COMPUTE  NEW  ANGLE  AND  ANGVEe  * 

AN«LE»OLDANG*ANGVEI.«OEL.TA  ( I )  *0.5*  A  NO  ACC*  (DELTA  I  I  )  **?) 

*  CHECK  FOR  NLG.  ANGLE, IE  SO  * 

PROCEED  TO  NEXT  DE°TH 

IF(ANGLE)  50*40,40 
CONT INUt 

call  veloc (phessr.watrho*ghavty*cosmax,anule» 

OEPTht I ) ,WEF0IA, AN3VEL) 

IF  (ANGLE. LE.OLi-ANG)  ANGVEL  —  ANGVEL 


A-B 


GO/1 


n  r  r  n  n  r  n  r,  r.  n 
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PROGRAM  CAVITY  73/74  OPT«l  FTN  4.6*45? 


OLOANG-aNGLE 


*  'fcOMFUTE  NEw  ANG.  ACC.  * 

CALL  ACCEL.  (PRES5H.WATRH0.GRAVTY.C0S  4AX  *  ANGLF  , 
.  DEPTH! I ) .REFDIA.ANGACC) 

ARTIME-AHTIME*OhLTA(I) 

IF ( ISTOwE-NSTORE( I) )  30.20.20 
■>u  CONTINUE 
STOP 


A***********************************. 

FORMAT  STATEMENTS 


lonu  format (bFio.si 
101U  FORMAT  (Fill. 5, 15) 

1020  FORMAT (2F10.5. 15) 

2000  FORMAT ( 1H1 .///////.45X. 

.  44H*  REFERENCE  INFORMATION*. 

.  //////, 47X.33HMQ0EL  CONFIGURATION  «  SUHHOC  .//. 

.  44X.2AHM0DEL  REF.  OIAM,  (IN.)  »  .F A . 2 . / / . 4HX , 

.  22HMOOEL  length  (in.)  *  .F5.2.//.5nrx. 

.  20HCAV.  DRAG  COEFF.  ■>  •  F5 . 3. // .  4  3*  . 

.  27HEFF  •  ORaG  AREA  (SO. IN.)  *  ,  FB  .5 , // 1 4BX ♦ 

.  22HM00EL  WtlGHT  <LflS>  ■  .FB. 5 . // .  40 X . 

.  30HMODEI.  DENSITY  ( LHS/CU  .  IN.  )  «  .FT  ,  4  . // .  4  1  X  , 

.  29HWATEH  OtNS,  ( SLUG/CU. FT . )  a  , FT . 5 , // , 42X  . 

.  28HATM.  PRESS.  (LOS/SQ.FT.)  a  . F 7  .  ?  , // . 4 1  X  , 

.  29HGHAV.  CONST.  (FT/SEC/SEC)  ■  , Fft . 3 » //  .  4 / X . 

.  23HENTRY  VtL.  (FT /SEC)  *  *F5 . 1 . // .56X . 

.  14HNO.  DEPTHS  *  .12*/) 

2005  FORMaTUH  .51  X,  IHhHETARD  (1/FT.)  ■  .  Fb .  5 .  /  ) 

2010  FORMAT  (1H  . 4M  X . 21 HDEPTHS  (CALIBERS)  a  .F4.1.4FA.1./) 

2 (lt'0  FORMAT  (1H1  .///////. 3SX, 

.b6H-  SPHERICAL  LAMINA  PERFORMANCE  DA 
.A  -.///, 53X.29HZEHU  TIME  REF.  *  WATEH  IMPACT) 

2010  FORMATUH  .////////. 34X.27H*  *************  , 

.  F4.1.38H  CAL  LAMINA  *************  #,///, 

.  5 1 » , 19HT IMt  STEP  (SEC)  *  *F6. 5 » / . AO  X , 1 OHNSTORE  *  , 

.  I3./.51X. 1 9dM AX  ANGLE  (DEG)  *  ,F5 ,?,///, 37X .4HA0S  .  , 

.  BX.THAPPROA.i 1SX.4HANG. .BX.4HASG. *  AX . 5HM00EL . / . 3 7 X  . 

.  4HTlME.BX*ttHCAV.DIA.i2X»SHANGl.  E»7X.4HVEL.»8X»4HACC«i 

.  BXi5HDEPTH»/«36X»5H(SEC>  » 0X , 7H (C AL»S ) . 3X ,5H (DEG l » 

.  4X.9H  iOEQ/SEC) ,3Xi 10H(DE3/SEC2) .4X,7H(CAL"S) ./) 

2040  FORMATUH  .3F.X.F5.4.  SX  .FS.2.4X  .F5 .2. 4X  *  F  A  ,  2, 1 X  ,F  1 4 . 1 . 3X  •  FT  .  2 ) 
END 


A-fl 


BO/  1 1 
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SuBHOUTl'tE  ACCEL 


73/74  opt  a  i 


FTN  4,h  +  45i> 


SUBROUTINE  ACCEL <  HHESSR  ,  W  ATRHO,*  GRAVTY*COSMAX*ANGLE* 

.  UtPTH,MEFDIA,ANGACC> 

COSAMO«COS(ANOUE) 

SINAm6«SI"U  ANGLE) 

A* (2.U-COSANO>#(l.U*COSAN3)/ (COSANG* < l.O-COSAMfl) ) 
HLOGX*A|_OG(X) 

A»2.i>-4.0*COSANG 

B»RLOOX*COsAMG* (<>.U*COSAN3**3-2.0*COSANG*COSANG-COSANU) 
Ca2.n*C0SA  JG/(SlNANG*5INANG) 

Z  = (A/ri) *C 

DUMMY*  (1,0*  (COSANu-COSMAX)  *Z)  /  <  SINANG**WLOf4K> 
ANGACC*-2.0*  (DEPTrt*RtFOIA*GRAm»PRESSR/rfATRHO>  *DUMMY/ 

.  ( <DEPTH*REFDIA>**2) 

heturn 

END 


SUBROUTINE  Vtuuc 


73/74  OPT  *  1 


FTN  4.b*452 


SUBROUTINE  VFLOC(PHESSR.WATRMO,QRAVTY*Cr)SMAX. ANGLE* 

.  DEPTH, REFDIA*ANGVEL) 

cosA‘iuacos  <  angle  ) 

SINAnG«SIN<  ANGLE) 

X*(2.U-C05rtNG>  *C1  .U  +  COSANG)  /  (COSANG*  (  1  .O-COSANG)  ) 

RLUGr ALOG ( X ) 

ANGVFL=SO’l (4.0* (utPTH»HEFOI A»GPAVTY+DRES?R/WATHHO) 

.  it  (COSANU-COSMAX)  /  (  (DEPTH*REFDTA) **?*SINANG 

.  <*S I NANO*RLOG )  ) 

RETURN 

END 


BO/K 


60/ 1  ( 
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•REFERENCE  INFORMATION 

MODEL  CONFIGURATION  a  RIGHT  CYL. 

MODEL  REF •  DIAM,  (IN-)  »  l.SO 
MODEL  LENGTH  (  IN.  )  ■  2,15 

CAV.  ORAG  COEFF.  «  .807 

tFF.  DRAG  AHEA  (SQ.IN.)  *  1.76715 

MODEL  WEIGHT  (L.BS)  *  1.07522 

MODtL  DENSITY  (LBS/CU.IN.)  «  ,2830 

WATER  DENS.  < SLUG/CU  .FT  . )  *  1.93946 

ATM.  PRESS.  (LBS/SU.FT . )  a  2116. HO 
GRAV.  CONST.  (FT/SEC/SEC)  *  32.174 

ENTRY  V EL.  (FT/SEC)  *  250.0 

NO.  DEPTHS  a  3 
RETARD  (1/FT.)  =  .28737 

DEPTHS  (CALIBERS!  «  10,0  20.0  30.0 


A-8 


NSWC  TR  81-69 

-SPHERICAL  LAMINA  PERFORMANCE  DATA 
ZERO  TIME  REF.  «  WATER  IMPACT 


•»««•*»•««•  *.*  10t0  CAL  LAMINA 


T V ME  STEP  (SEC)  =  .00010 

NSTORF  *  5 

MAX  ANGLE  (DEG)  >=  9,45 


ARS. 

APPROX. 

ANfi. 

ANG  . 

MODEL 

T 1  ME 

C  A  V  ,  0 1  A  , 

ANGLE 

VFL  . 

ACC. 

DEPTH 

<SEC) 

<CAL"S) 

(DEG) 

(DEG/SEC) 

(DE6/SEC2) 

(CAL“S) 

.nO60 

1.00 

2.86 

2524.23 

-2148373,6 

10.00 

.0066 

1, 37 

3.93 

1837.22 

-911064.9 

10.69 

.0070 

1.66 

4.75 

1486.39 

-548801,7 

11.36 

.0  0-75 

1.89 

5.43 

1257.18 

-386432.5 

12.02 

.0080 

2.10 

6.02 

10RR.26 

-297256.2 

12.66 

.0085 

2.27 

6.53 

954.43 

-242139.3 

13.29 

.0090 

2.43 

6.97 

843.15 

-205324.0 

13.90 

.0095 

2.57 

7.37 

747.34 

-179366.2 

14.50 

.0100 

2 .69 

7,72 

662.66 

-160333.9 

15.09 

.0105 

2.80 

8.03 

586.25 

-145974.0 

15.67 

.0  110 

2.89 

8.31 

516.15 

-134912.2 

16.  ?3 

.0115 

2.97 

8.55 

450.94 

-126269.1 

16.78 

.0120 

3.05 

8. 76 

369.58 

-119459.8 

17.32 

.0  125 

3.11 

8.94 

331.25 

-114003.9 

17.86 

.0130 

3.16 

9,09 

275.31 

-109661.4 

18.38 

.0135  ' 

3.20 

9.22 

221 .24 

-106594.0 

18.89 

.01  AO 

3.24 

9.31 

168.59 

-104140.8 

19.39 

.0  145 

3.26 

9.39 

1 16.9b 

-102403.2 

19.89 

.0150 

3.28 

9,43 

66.08 

-101314.7 

20.37 

.0155 

3.28 

9.45 

15.57 

-100835.6 

20.85 

.0160 

3.26 

9.45 

-26.66 

-100883.8 

21.32 

.0  165 

3.28 

9 1  42 

-76.22 

-101462.8 

21.78 

.0170 

3.26 

9.37 

-127.24 

-102697.1 

32.23 

.0175 

3.23 

9.30 

-179.02 

-104571.7 

22.68 

.0180 

3.20 

9.19 

-231.93 

-107178.7 

23.12 

.0  185 

3,15 

9.07 

-2B6.34 

-110624.3 

23.55 

.0190 

3.10 

8.91 

-342.71 

-115059.4 

23.98 

.0195 

3.03 

8.72 

-401.59 

-120696.5 

24.40 

.0200 

2.9b 

8.51 

-463.65 

-127835.9 

24.81 

.0205 

2.87 

8.26 

-529.74 

-136908.7 

25.22 

.0210 

2.77 

7.96 

-bOO .98 

-148547.8 

25.62 

.0215 

2.66 

7.66 

-678.86 

-163711.7 

26.01 

.0220 

2.54 

7.30 

-765.51 

-183911.0 

26.40 

.o2?5 

2.40 

6.69 

-864.00 

-211650.  6 

26.  79 

.0230. 

2.24 

6.43 

-979.11 

-251366.2 

27. IT 

.0  235 

2.06 

5.91 

-1118,70 

-311619.6 

27.54 

.n240 

1.85 

5.30 

-1297.04 

-411041.0 

27.91 

.0245 

1  .60 

4.60 

-1543,69 

-597968.7 

28.28 

.0250 

1.30 

3.74 

-1934.54 

-1040086.3 

28.64 

.0255 

.91 

2.60 

-2769.24 

-2802135.3 

28.99 

.0260 

.12 

.34 

-862860962.9 

29.34 
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*  •  4  *  # 

#  #  *  »  * 

*»*  20.0  CAL  LAMINA 

•o#  «•**** 

»  *  #  #  i 

TIME 

STEP  (SEC)  ■ 

.00010 

NSTORE  ■ 

5 

MAX  ANGLE  (DEG)  ■ 

3.24 

ARS. 

APPROX. 

ANG. 

ANG. 

MODEL 

TIME 

CAV.OIA 

.  ANGLE 

VfL. 

ACC. 

depth 

(SEC) 

( C  AL'*S ) 

(OEG ) 

(DEGySEC) 

(0EG/SEC2) 

<CAL"S> 

.0146 

1.00 

1.43 

760,50 

-455800.2 

20.00 

.nisi 

1.23 

1.77 

590,90 

-257576.3 

20.48 

.0156 

1.42 

2.03 

484.92 

-176413.0 

20.96 

.0161 

1.S7 

2.26 

408.34 

-133996.1 

21.43 

.0166 

1.71 

2.44 

348.19 

-108615.4 

21. 8t) 

.oin 

1,82 

2  a  60 

298.29 

-92084.1 

22*34 

.0176 

1.91 

2.74 

255 .27 

-80692.0 

22.78 

.0101 

2.00 

2.86 

217.08 

-72541.1 

23.22 

.0186 

2.07 

2.96 

18?. 38 

“■66566  •  6 

23.65 

.0191 

2.12 

3.0  4 

150.26 

-62136.2 

24.07 

.0196 

2.17 

3.11 

120.06 

-58854.9 

24.49 

.  0201 

•  d  1 

3.16 

91  .26 

-56469.1 

24.90 

.0206 

2.23 

3.20 

63.47 

-54814.1 

25.31 

.0211 

2.25 

3.23 

36.35 

-53785.4 

25.71 

.0216 

2,26 

3.24 

9.60 

-53321.9 

26.10 

.0221 

2.26 

3.24 

-14.93 

-53371.1 

26.49 

.0226 

2  ■  25 

3.22 

-41,73 

-53944.5 

26.80 

.0231 

2.23 

3.20 

-68.97 

-55092.4 

27.25 

•  0236 

2.20 

3.15 

-96.93 

-56883.6 

27.63 

.0241 

2.16 

3.10 

~i25.97 

-59433.1 

28.00 

.0246 

2.11 

3.03 

-156,51 

-62920.4 

28.36 

•  0  251 

2.05 

2.94 

-189,09 

-67622.3 

28.72 

.0  256 

1.98 

2.64 

-224.40 

-73971.5 

29.07 

.0261 

1.90 

2.72 

-263.44 

-B2667.5 

29.4? 

.0266 

1.80 

2.57 

-307.64 

-94896.4 

29.77 

.0271 

1.68 

2.41 

-359.26 

-U2B15.1 

30.11 

.0276 

1,54 

2.21 

-422.07 

-140719.8 

30.45 

.0281 

1.38 

1.98 

-503.16 

-188392.3 

30.70 

.0286 

1.19 

1.71 

-617.90 

-283058.5 

31.11 

.0291 

.94 

1.35 

-809.73 

-532201.3 

31.43 

.0296 

.60 

.85 

-1298,35 

-1912956.5 

31.75 
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30.0  CAL  LAMINA 


TIME  STEP  (SEC  1  » 
NSTORE  ■ 
MAX  ANGLE  (DEG)  ■ 

.00010 

5 

1.48 

ARS  . 

T  IMF 
<  SEC ) 

APPROX. 

CAW.OIA.  ANGLE 
<CAL,(S)  <  OEG ) 

ANG  . 

VFL. 

(0EG/5EC) 

ANG. 

ACC. 

t  DEG/SEC2) 

MODEL 

DEPTH 

<CAL"S) 

0270 

1.00 

0275 

1.13 

0280 

1.24 

02R5 

1.33 

0290 

1.40 

0295 

1.46 

0300 

1.50 

,0305 

1.53 

.0310 

1.55 

.0315 

1.55 

.0320 

1.55 

.0325 

1.53 

.0330 

1  .50 

,0335 

1.46 

.0340 

1.40 

.0345 

1.33 

.0350 

1.25 

.0355 

1.14 

.0360 

1.01 

.0  365 

.63 

.0370 

.60 

.0375 

.15 

.95 

2BR.B2 

1.06 

230.10 

1.19 

166.17 

1.27 

150.56 

1.34 

120.12 

1.39 

93.05 

1.43 

68.22 

1.46 

44.87 

1.4B 

22.42 

1.48 

.41 

1.48 

-21.61 

1.46 

-44.03 

1.43 

-67.34 

1.39 

-92.10 

1.34 

-119.08 

1.27 

-149.37 

1.19 

-184.75 

1.09 

-22P.29 

.96 

-286.29 

.80 

-374.45 

.57 

-552.06 

.14 

-2170.53 

-140130. A  30.00 

-99634.3  30.34 

-78150.5  20.67 

-65303.5  31.00 

-57071,8  31.33 

-61608.7  31.65 

-47969.1  31.97 

-45637.5  32.29 

-44328.6  32.60 

-43697.3  32.91 

-44297.9  33.21 

-45572.7  33.51 

-47863.0  33.81 

-61447,5  34.11 

-56830.8  34.40 

-64936.6  34.69 

-77564.2  34.97 

-98613.0  35.26 

-138067.7  35.54 

-229697  *  3  35.82 

-575366.9  36.09 

-31522971.7  36.36 
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APPENDIX  D 


nPRIVATION  OF  LAMINA  MOTION  EQUATIONS  WHICH  INCLUDE 
derivation  Ot  n„PHKRTC  SURFACE  pressure 


LINE  OF  FIRE 


NSWC  TR  81-59 


This  derivation  uaee  an  energy-approach  to  obtain  the  equations  governing 
the  motion  of  the  lamina.  Refer  to  the  figure  above  throughout  the  derivation. 

First,  obtain  expression  for  the  potential  due  to  surface  pressure  P0, 
which  is  equal  to  negative  of  the  work  done  by  P0  moving  thru  the  displacement 
(Qb  -  90°) : 


.  dV„  -  -  dWD  -  -  (-  P0  A,  r  d9b)  (1) 

ro  ro 

where  the  area  Aa ,  as  shown  in  the  figure,  is  given  by: 

Ag  -  2%  Ar  sin  8b 


Substitute  expression  for  Afl  into  Equation  (1)  and  integrate  to  obtain  an 
equation  for  Vp  in  terms  of  8b,  thus: 


PQ  21T  r2  Ar  sin  8b  d0b 
90° 


2tt  r2  Ar  p0  cos  0  b 


90° 


or: 


Vp  *  -  2TT  r2  Ar  P0  cob  0b  (2) 

o 

It  can  be  shown  that: 

cos  0  b  ■  cos  0  a  -  1 


Subatituting  this  into  Equation  (2)  yields  the  expression  for  Vp  in  terms  of 
angle  0a:  ’  c 


Vp  -  -  2  r2  r  P0  (cos  a  -  1)  (3) 

Next  consider  gravity  potential  Vg.  If  the  original  undisturbed  water  surface 
is  taken  as  the  reference  datum  for  zero  gravity  potential  then 
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Vg  ■  -  2n  r2  Ar  p„  g  h 
where,  as  shown  in  figure: 

h  ■  r ( 2  cos  9a  -  1 ) / 2 

Hence  gravity  potential  in  terms  of  shell  angle  is 
Vg  ■  TTr-^  Pw  Ar  g  (1  -  2  cos  Ga) 

Finally,  combine  potentials  due  to  gravity  and  surface  pressure  to  obtain  the 
total  potential  of  the  shell/pressure  system 


V  ‘  vg  ♦ 


3 

TTr  p  Ar  g  (1  -  2  cos  0  )  -  27T  r2  Ar  P  (cos  0  -  1) 

V  A  O  A 


(4) 


Next,  consider  the  kinetic  energy  of  the  moving  shell, 
shell  angle  0a,  the  expression  for  the  kinetic  energy  can  be 


In  terms  of  the 
shown  to  beA-l 


Tf 


y  r^  Ar  8a  sin2  8a  In 


(2  -  cos 


V 


(1  ■»  coa 


V 


cos  0.  (1  -  COB 


V 


(5) 


For  the  total  energy  in  the  conservative  shel l/surface-preseure  system  the 
following  relation  holds  at  all  times 


(Tf  +  V)f  -  (Tf  +  V)2  -  Const. 


(6) 


When  the  shell  is  at  the  peak  of  its  swing,  its  kinetic  energy  is  zero  and  the 
total  potential  of  the  system  is  at  a  maximum  which,  from  Equation  (6),  is  equal 
to  the  total  energy  of  the  shell  at  any  other  time. 

(Tf  +  V)  -  V(peak)  .  (7) 


At  the  peak  position  9a  is  also  at  a  maximum  and  is  designated  as  0ms 
the  maximum  value  that  9fl  ever  reaches.  When  the  previously  derived 
expressions  for  Tf  and  V  are  substituted  into  Equation  (7)  and  the  shell  angle 
in  the  expression  for  V(peak)  is  replaced  by  9m,  Equation  (7)  becomes 


A-lAbelson,  H. ,  "Cavity  Shapes  at  Vertical  Water  Entry  -  A  Comparison  of 
Calculated  and  Observed  Shapes,"  NOLTR  67-31  (1967). 
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^  r  ^  (Ar)  6*  a in^  0a  In  [x]  +  it  (Ar)  ft,  g  (1  -  2  co*  0^) 
~  2 it  r2  (tr)  ?0  (cos  0a  -  1) 

■  H  pw  (Ar)  g  (1  -  2  coa  0m)  -  2tt  (Ar)  PD  (coa  0m  -  1) 


where 


X  »  (2  -  cos  0a)  (1  +  coa  0 a ) / coa  0a  (1  -  coa  0a) 


Solving  this  last  expression  for  the  angular  veloctiy  9 a  yields  the  governing 
differential  equation  for  the  shell  angle 


ft 

rg  +  (p0/Pw5 

0 

COl  0  -  COB  0 

a  m 

»in2  0  In 

s 

’(2  -  coa  0a)  (l  +  cob  0*)"] 

coa  0a  (l  -  coa  0a)  j 

(8) 


The  expression  for  angular  acceleration  9a  can  be  found  by  differentiating 
Equation  (8)  with  respect  to  time 

-2  Jrg  +  (PQ/PW^|  ["  1  +  (cob  0a  -  cos  0m)  CZ) “| 

e*  “  J  ? —  L  nrcT"  w — ~J 


(2  -  coa  6a)  (1  +  cos  9a) 
cos  9a  (1  -  cos  9a ) 


(9) 


Z  ■ 


2  -  A  cos  9a 


2  cos  Ba 


cos  0a  (cos 3  0^  r'2  cos2  8a  -  cos  9a  +  2)  In  [xj  sin2  0a 


The  expression  for  cos  9m  is  derived  next.  At  the  peak  of  its  outward  swing 
(corresponding  to  0a  “  9m)  all  the  kinetic  energy  imparted  initially  to  the 
lamina  by  the  passing  projectile  has  been  converted  entirely  into  potential 
energy  thus: 

vpeak  "  ^iri  “  ^r^ 

where  Tp,  the  kinetic  energy  of  the  projectile,  is  given  by 


1  2 
TP  "  2  nlP  VP 
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If  Che  cavity  drag  coefficient  is  assumed  to  be  constant  then  Tp  can  be 
written  in  terms  of  lamina  depth  r  as 


i  mp  vp2  e 
1  o 


-2ar 


The  expressions  for  Tp,  Vpea|<  and  V^n  are  substituted  into  Equation  (10) 
and  the  differentiation  carried  out  to  obtain 


IT  P  p  (£r)  g  (1  -  2  c os  0„,)  -  (-it  r-3  p  (Ar)  g) 


w 


“  2 it  rZ  (Ar)  P0  (cos  Q_  -  1)  -  a  m„  vD  e  2Qlr(Ar) 


P  yPr 


Solving  for  cos  6-,  yields  the  final  desired  relation 


cos  9m  -  j 


a  m  vn  e 
p  Ho 


-2ar 


"  r2  <r  Pw  8  +  P0>j 
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